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Abstract This article reviews recent studies of scale interactions in magnetohydrodynamic 
turbulence. The present day increase of computing power, which allows for the exploration 
of different configurations of turbulence in conducting flows, and the development of shell-to- 
shell transfer functions, has led to detailed studies of interactions between the velocity and the 
magnetic field and between scales. In particular, processes such as induction and dynamo action, 
the damping of velocity fluctuations by the Lorentz force, or the development of anisotropics, 
can be characterized at different scales. In this context we consider three different configurations 
often studied in the literature: mechanically forced turbulence, freely decaying turbulence, and 
turbulence in the presence of a uniform magnetic field. Each configuration is of interest for 
different geophysical and astrophysical applications. Local and non-local transfers are discussed 
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for each case. While the transfer between scales of solely kinetic or solely magnetic energy 
is local, transfers between kinetic and magnetic fields are observed to be local or non-local 
depending on the configuration. Scale interactions in the cascade of magnetic helicity are also 
reviewed. Based on the results, the validity of several usual assumptions in hydrodynamic 
turbulence, such as isotropy of the small scales or universality, is discussed. 
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1 Introduction 

Turbulence is a multiscale phenomenon ubiquitous in geophysical and astrophys- 
ical flows. In many of these flows, the coupling of a conducting fluid with elec- 
tromagnetic fields requires consideration of the magnetohydrodynamic (MHD) 
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equations (see, e.g. 



T 

Moffatt 197a ) . The equations describe the dynamics of non- 



relativistic conducting fluids as, e.g., in the Earth's core or in industrial appli- 
cations, and under some approximations they can also describe the large-scale 
behavior of magnetospheric, space, and astrophysical plasmas. In these latter 
cases, care must be taken to consider only the scales where a one-fluid approxi- 
mation holds, as scales small enough may require consideration of kinetic plasma 
effects such as ambipolar diffusion in weakly ionized plasmas as the interstellar 
medium, or the Hall current for highly ionized media such as small scales in the 
solar wind. However, in those cases the MHD equations still give a good de- 
scription of large scales, and the approximation gives a useful approach to get 
lowest-order physical insight into the fate of the flows. 

In the simplest case, that of an incompressible flow with constant mass density, 
the equations give the evolution of the bulk fluid velocity u and of the magnetic 
field b: 

d t u + u • Vu = - Vp + b • Vb + zA7 2 u, (1) 

d t b + u • Vb = b • Vu + r/V 2 b, (2) 

where the magnetic field is written in Alfvenic units, the density is set to unity, 
and p is the (fluid plus magnetic) pressure. The kinematic viscosity v and mag- 
netic diffusivity r] control respectively the viscous and Ohmic dissipation. These 
equations are constrained by the incompressibility condition and by the solenoidal 
character of the magnetic field, 

V-u = 0, V-b = 0. (3) 

Two different Reynolds numbers can be defined in MHD flows. The mechanical 
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Reynolds number 

Re=^, (4) 

which is the ratio of convective to viscous forces (where U is the rms velocity and 
L a characteristic lengthscale of the flow), and the magnetic Reynolds number 

Rm = — (5) 
rj 

that can be interpreted as the ratio of induction to Ohmic dissipation. In many 
flows these Reynolds numbers are very large, and the flows are in a turbulent 
regime. 

While in the hydrodynamic case the phenomenological theory of Kolmogorov 
(K41) gives to a good approximation (albeit intermittency corrections) the power 
law of the energy spectrum, no clearly established phenomenological counterpart 
exists in MHD. This has many implications as the energy dissipation rate (re- 
quired to predict, e.g., heating rates in solar and space physics) depends on the 
slope of the energy spectrum. Also, subgrid models, required to do numerical 
modeling in astrophysics and geophysics given the large scale separation involved 
in such flows, are less developed in MHD as a result of the lack of detailed knowl- 
edge of its energy spectrum. 

In the Kolmogorov description of hydrodynamic turbulence, the interactions of 
similar size eddies play the basic role of cascading the energy to smaller scales on 
a scale-dependent time scale Tg ~ £/ug, where t is the examined length scale and 
ug the characteristic velocity at that scale. This time scale, which is proportional 
to the eddy turnover time at the scale £, is the only time scale available on di- 
mensional grounds in the inertial range, provided enough scale separation exists 
between forcing and dissipation. In this context, interactions between scales are 
local (in spectral space) as dominant interactions are between eddies of similar 
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sizes. One then expects the statistical properties of sufficiently small scales to 

be independent of the way turbulence is generated, and to have therefore uni- 
versal character. Recent experiments showed deviations from this behavior even 



or presence of long-time correlations in the small scales; see 


Car 


lier et al. 2001. 


Poulain et al. 2006 


Shen & Warhaft 200ol. Iwiltse & Glezer 1993. 


1998 


). Numer- 



ical simulations also gave evidence of the presence of non- local interactions wit l 
the large scale flow playing a role in the cascade of energy (lAlexakis et al. 2005al . 



Domaradzki 1988 



Domaradzki k, Rogallo 1990i . 



Zhou 19931 ). In numerical sim- 



ulations with Reynolds numbers as high as R\ ~ 800, it was observed that 
20% of th e energy flux in the small scales was due to interactions with the large 



scale flow (jMininni et al. 20061 ). However, more recent simulations with Reynolds 
numbers up to R\ ~ 1300 using spatial resolutions of 2048 3 grid points showed 
that as the Reynolds number is increased, the percentage of the non-local flux 
decreases as a power law of the Reynolds number, suggesting that the flux in hy- 
drod ynamic turbulence m ay be predominantly local for very large Reynolds num- 



ber 



Mini nni et al. 20081). Recent theoretic al results put this in firmer grounds 



(lAluie & Evink 2009al . 



Evink &: Aluie 20091 ) . showing that the energy flux in hy- 



drodynamic turbulence is local in the limit of infinite Reynolds number and 
obtaining bounds on the scaling of the non-local contribution to the flux with 
Reynolds number which are in agreement with the numerical results. 

The case for MHD turbulence is less clear and has given rise to more controver- 
sies. Several attempts have been done to ext end the phenomenological arguments 



of Kolmogorov to conducting flow s (see, e.g..lBoldvrev 2006 



Iroshnikov 1 963. 



Kraichnan 1965 



Matthaeus fc Zhou 19891 ). However, the MHD 



iGoldreich fc Sridhar 1995 
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equivalent of the 4/5 law in hydrodyna mic turbulence (the Politano-Pouquet re- 



lations, see 



Politano Pouquet 1998a 



bj) couple the velocity and the magnetic 



field in a way that can be compatible with several power law behaviors; in three 
dimensions they read 

± , n ,2\ 4 



5zf(l) 1^(1)1 



(6) 



where are the dissipation rates of the Elsasser variables z ± = u ± b, and the 
sub index || denotes the increment of the field along the displacement vector 1. 

Moreover, even in the simplest incompressible case, at least two time scales 
can be identified in the inertial range of MHD turbulence. Besides the eddy 
turnover time, incompressible MHD flows are also characterized by the period 
of Alfven waves r ~ (BqL) -1 , where Bq is the amplitude of the large scale 
magnetic field in Alfvenic u nits. In a first attemp t to derive a phenomenologi- 



cal theory, llroshnikov (19631 ) and 



Kraichnan (19 65) (IK) assumed that the large 



scale magnetic field acts as a uniform field for the small scale fluctuations, which 
then behave as Alfven waves. In that case, small scales can interact not only 
through the eddies but also through Alfven packages, which reduce the energy 
flux to small scales by increasing its transfer time. This introduces in prac- 
tice a non-local interac tion as the waves propagate along the large scale field 



sec 



Gomez et al. 19991 for a discussion). From dimensional analysis, Iroshnikov 
and Kraichnan then derived an isotropic energy spectrum proportional to £;~ 3 / 2 . 
Later, extensions where considered to take int o account the anisotropy induced at 



sma 



2005 



1 scales by the large scale m agnetic field (jBoldvrev 20061 . 



Galtier et al. 2000l . 



Goldreich fc Sridhar 19951 ). Some of these extensions, after accounting for 
the anisotropy, rely on some form of a balance between the two fields that leaves 
only the turnover time as the relevant time scale, and can be therefore considered 
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local or non-local depending on the authors. 

At the core of the early disquisitions is the fact that in MHD the roles of a 
large scale flow and of a large scale magnetic field are different. While a (uni- 
form) large scale flow can be removed by a Galilean transformation, a large scale 
magnetic field cannot. As a remarkable coincidence, lack of Galilean invariance 
is at the basis of the ~ k~ 3 ^ 2 spectrum for hydrodynamic tur bulence within th e 
framework of the Direct Interaction Approximation (DIA) bv lKraichnan (19591 ). 
a flaw later corrected by the development of the Test Field Model (TFM) and 
the Lagrangian History DIA (LHDIA). However, in MHD magnetic fields are not 
Galilean invariant and for this reason the associated Alfven waves have to be taken 
into account in phenomenological theories, and are also considered when studying 
non-l ocal effects in the E ddy-Damped Quasi-Normal M arkovian (EDQNM) clo - 



sure (jPouquet et al. 19761 ) or in weak turbulence theory (jNazarenko et al. 2001 
However, although phenomenological descriptions assume that a large scale field 
has the effect of reducing the energy cascade rate, the transfer of energy (and the 
cascade) in many cases still takes place between eddies of similar size, presumably 
allowing for recovery of universal statistical properties at small scales. 

In recent years, this universal behavior has been questioned by different au- 
thors. Because energy can be injected in MHD by a mechanical forcing or by 
an electromotive forcing, MHD turbulence is characterized by a larger number 
of regimes than hydrodynamic turbulence even in its simplest configurations. 
Magnetic fields in planets and stars are believed to be generated by dynamo 
actio n, where turbulent motions sustain magnetic fields against Qhmic dissipa- 



tion ( Brandenb urg Subramanian 2 005 



Krause &: Raedler 198C 



Moffatt 197S 



Pouquet et al. 19761 ). This regime is often studied in numerical simulations (and 



recently in experiments, see 



Mininni 

Monchaux et al. 20071 ) by mechanically stirring the 



flow. Depending on the amount of mechanical helicity in the flow (the alignment 
between the velocity and the vorticity), or in the presence of large-scale shear, 
the magnetic field generated may have large or small-scale correlation (compared 
with the integral scale of the flow), giving a steady state that may be dominated 
by mechanical or magnetic energy. On the other hand, plasmas in the solar corona 
and in the solar wind are dominated by magnetic energy, and are often studied 
numerically by stirring the flow with electromotive forces or using simulations 
of freely decaying turbulence. Finally, the amount of cross-correlation between 
the velocity and the magnetic field depends on the flow (e.g., on the heliocentric 
distance in the solar wind) and can also be varied in the simulations. 

The questioning of universality was accompanied by recent detailed studies of 
scale interactions in MHD turbulence. Many of the studies considered the so 
called shell-to-shell transfer functions and partial energy fluxes, either in numer- 
ical simulations, observations, closures, or from the theoretical point of view. In 
the following sections we give a review of the results in this area, considering the 
several regimes studied by different authors, and also some examples of possible 
sources of non-locality in MHD. Finally, we discuss the results in the context 
of universality and of phenomenological theories for MHD. To briefly summarize 
the results, several authors have shown that the locality of energy transfer is in 
question in MHD flows. In particular, it was shown from simulations that the 
transfer of energy in MHD has two components: a local one that shares similar 
properties with hydrodynamic turbulence, and a component coupling the veloc- 
ity and magnetic fields for which energy from the large scales can be under some 
circumstances injected directly into the small scales without the intervention of 
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intermediate scales. 
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2 Indirect evidence of non-locality 

Some theoretical, phenomenological, and (more recently) numerical results in- 
dicate scale interactions in MHD can be, under some conditions, of a different 
nature than in hydrodynamic turbulence. In this section we review early theoret- 
ical indications of non-locality in MHD turbulence, as well as numerical results 
that support the theoretical arguments without directly measuring scale interac- 
tions. 

Early studies of dynamo action, and of magnetic field evolution under flows 
with simple strain, show that a large scale flow can excite, through field line 
stretching, magnetic fields at widely separat ed scales. One of the first works along 



these lines is the work of 



Batchelor (19501 ) where he considered the similarity 



between the induction equation and the vorticity equation (u> = V x u): 

— + u • Vu = u) ■ Vu + uV z oj. (7) 
at 

While the second term on the l.h.s. advects the vorticity, the first term on 
the right (in three dimensions) produces vorticity by vortex stretching. For 
Pm = n/u > 1 (the magnetic Prandtl number) Batchelor then concluded that 
the magnetic field would grow as magnetic field line stretching overcomes Ohmic 
dissipation. Later works considering stretchin g by uniform straining motion 



dMoffatt k Saffman 1964 



Zel'Dovich et al. 1984 ) showed that a lar ge scale mag- 



netic field can directly create small-scale magnetic fields. The work of lKazanstev (1 968) 
considered a similar process under a random velocity field and described a non- 
local coupling that sustains the so-called small-scale dynamo, where magnetic 
fields are amplified at scales smaller than the integral scale of the flow. Several 
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numerical simulations support these results, and show that smooth motions at 



the viscous scale give ex ponential growth of magne t ic fie 



magnetic diffusion scale (jSchekochihin et al. 2002al . 



2004 



ds that can peak at the 



2002bh . 



The opposite limit, when the magnetic Prandtl number is much smaller than 
unity (a case of interest fo r industrial flows), is some times studied using the quasi- 



static approximation (see 



Knaepen fc Moreau 20081 for a review). In this case, an 



external uniform magnetic field is applied, and the magnetic Reynolds number 
is chosen small enough that magnetic field fluctuations are rapidly damped. In 
that limit the Lorentz force in the momentum equation reduces to linear Joule 
damping 



<9u „ 



2V7-2 



-Vp + aB^V 



u 



+ vV 2 u, 



(8) 



where a is the conductivity of the medium and the uniform magnetic field B$ is 
in the z direction. The Joule damping, although anisotropic in spectral space, is 
roughly independent of the wavenumber, and unlike viscous damping is not con- 
centrated at small scales but rather acts at all scales. As a result, the large-scale 
magnetic field in this approximation exerts work over all scales in the velocity 
field (damping turbulent fluctuations) in a non-local way. We will see that the 
shell-to-shell transfers indicate in some cases similar behavior of the Lorentz force 
even in cases far from this approximation. 

Another important example concerns Alfven waves, which are also non-linear 
solutions of the ideal MHD equations. Alfvenic states with u = ±b make the 
non-linear terms in Eqs. ([T|) and (|2|) zero, leaving only interactions with the large 
scale fields to transport energy across scales. Finally, it is worth mentioning 

here some recent attempts to build shell mode l s of MHD turbulence (see e.g. 
1 p 1 



Lessinnes et al. 20091 . iPlunian &: Stepanov 20071 . 



Stepanov &; Plunian 20081 ). In 
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these models, it was found that many features of steady state MHD turbulence 

can be reproduced using local coupling between shells, but that to reproduce the 

small-scale dynamo and turbulence a t Pm ^> 1, non-local transfers have to be 



considered (jStepanov Plunian 20081 ). 



3 The shell-to-shell transfer 



In recent years, the increase in computing power allowed numerical exploration 
of MHD t urbulence in differen t regimes. The development of shell-to -shell trans- 



fers (see 



Alexakis et al. 2005b 



Par et al. 2001 



Debliquv et al. 20 05) allowed for 



explicit computation of detailed scale interactions in MHD turbulence using the 
output stemming from the simulations and without the need to compute the more 
expensive triadic interactions. In this section we briefly introduce the isotropic 
shell-to-shell energy transfer functions, and describe how fluxes can be obtained 
from them. 



A shell filter decomposition of the two fields is introduced as 



u ( x ) = H u ^( x )> 

K 



b(x) = £Mx), 



A' 



where 



and 



Ujf (x) = £ G ( k ) e4k ' X 

Ai<|k|<A 2 



Mx) = Y, b(k)e ik ' x 

Ai<|k|<A 2 



(9) 
(10) 

(11) 
(12) 



Here u(k) and b(k) are respectively the Fourier transforms of the velocity and 
magnetic fields with wavenumber k. The shell filtered fields and hx are 
therefore defined as the field components whose Fourier transforms contain only 
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wavenumbers in a given shell K. These shells can be defined with linear binning 

using K\ = K and K2 = K + 1, or alternatively with logarithmic binning using 
K\ = 7 n -K"o and K2 = 7 n+1 i£o for some positive 7 > 1 and integer n (7 = 2 is 
often used) . The latter definition has the advantage of being conceptually closer 
to the idea of "scale" of eddies in turbulence, which in general implies logarith- 
mic division of wavenumbers. The former has the advantage of having a direct 
association with Alfven waves, which are of the form u = ±b ~ e »(k-x±wt) j n 
periodic boxes or in infinite domains, and which are more akin to linear treat- 
ment of spectral space. Note that the transfer among logarithmic shells can be 
reconstructed by summing over the linearly spaced shells. 

Another variant when defining the shell filter decomposition has to do with 
the choice of usin g sharp filters (as in the equations above) or smooth filters 



(|Evink 1994 



2005). This issue has raised some controversy in the hydrodynamic 



case, with claims that non-localities observed in simulation s may be due to the 



commonly used sharp filters. Recent numerical comparisons (jDomaradzki fc Carati 2007a 
have shown that results are only weakly dependent on the shape of the filter used, 
except in the case where a very broad smooth filter is considered. Moreover, re- 
cent theoretical results were able to show locality of hydrodynamic turbulence 
in Fourier sp ace in the limit of infinite Reynolds numb ers for both smooth and 



sharp filters (|Aluie & Evink 2009a 



Evink & Aluie 20091 ) . 



Based on the shell filter decomposition, the evolution of the kinetic energy in 
a shell K, E U (K) = J u^/2cfa; 3 , can be derived from Eq. (P) as 
dE u {K) 



dt Q 



[Tuu(Q, K) + T bu (Q, K)\ - vD u (K), (13) 
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and for the magnetic energy, E b (K) = \b\j2dx 3 , from Eq. ([2]) 

^9 = [T bb (Q, K) + T ub (Q, K)\ - rjD b (K), (14) 
Q 

where the functions D U (K) and D b (K) express respectively the kinetic and 
magnetic energy dissipation in the shell K. The transfer functions T UU (Q,K), 
T ub (Q, K), T bb (Q,K), and T bu (Q,K), that express the energy transfer between 
different fields and shells are given by 

T UU (Q,K) = - J u K (u-V)u Q dx 3 , (15) 

T bu (Q,K) = Ju K (h-V)b Q dx 3 , (16) 

T bb (Q,K) = - Jb K (u-V)b Q dx 3 , (17) 

T ub (Q,K) = Jb K (h-X7)u Q dx 3 . (18) 

The function T UU (Q, K) measures the transfer rate of kinetic energy in the shell 
Q to kinetic energy in the shell K due to the advection term in the momentum 
equation flT]). This is the non-linear transfer that is also present in hydrodynamic 
turbulence. Similarly, T bb (Q, K) expresses the rate of magnetic energy transferred 
from the shell Q to magnetic energy in the shell K due to the magnetic advection 
term. The Lorentz force is responsible for the transfer of energy from the magnetic 
field in the shell Q to the velocity field in the shell K, as measured by T bu (Q, K). 
Finally, the term responsible for the stretching of magnetic field lines, the first 
term on the r.h.s. of Eq. (|2J), results in the transfer of kinetic energy from shell 
Q to magnetic energy in shell K and is expressed by T ub (Q, K). This is the term 
that describes magnetic induction and dynamo action. 
All these transfer functions satisfy 

T VW (Q,K) = -T WV (K,Q). (19) 
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(where the subindices v and w stand for u or b) . This expression indicates that the 

rate at which the shell Q gives energy to the shell K is equal to the rate the shell 
K receives energy from the shell Q, and is a necessary condition to define shell- 
to-shell transfers that satisfy a detailed energy balance between shells. Then, the 
contribution of these transfers to the total energy flux can be computed as: 



U mo (k) = -Y l ^2 T vw(K,Q). 

K=0 Q 



(20) 



Besides the total energy, the MHD equations have two more ideal invariants: 
the cross- helicity C = J u • b dx 3 , and the magnetic helicity H = J a • b dx s where 
a is the vector potential such as b = V x a. These quantities also satisfy detailed 
balance equations equivalent to Eqs. (|13p and fl!4|). Shell-to-shell transfer func- 



tions for the magnetic helicity have been defined in 
transfer from shell Q to shell K is measured by 



Alexakis et al. (2009 ). Its 



T H (K, Q) = Jh K - (u K x bg) dx 3 



(21) 



The energy transfer functions were also generaliz ed in recent works to consider 
the flux of energy in terms of the Elsasser variables (lAlexakis et al. 2007al.l2 005b 



Carati et al. 20061 ). anisotropic transfers (jAlexakis et al. 2007al . 



Teaca et al. 2003 ). 



forward and backwa rd transfers in an attempt to quantify 



for subgrid models (jCarati et ak 2006. 



jackscatter required 



Debliauv et al. 20051 ). and extensions to 



consider compressibility effects iGraham et al. f 20101) . and kinetic plasma effects 



as in two-fluid MHD approximations (jMininni et al. 2 007). 



4 Direct studies of multi-scale interactions 



The shell-to-shell energy transfers have been studied extensively (lA lexakis et a 



Carati et al. 2006 



Par et al. 2001 



Debliauv et al. 2 005 



Mininni et al. 2005al , 



. 2005b, 



Verma 20041 ) 
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for a variety of mechanically forced and decaying MHD flows in two and three 

dimensions. Depending on the configuration, different degrees of non-locality 

were reported. In the following subsections we present a short summary of 

the results discriminating by the forcing configuration. Overall, we can say 

that in all cases examined in the literature the transfers T uu and T^b have a 

local behavior: energy is transferred forward between nea rby shells, in a fashion 



similar to what is observed i n hydrodynamic turbulence (Alcxak is et al. 2005a 



Domaradz 



ri Rogallo 19901 . 



Mininni et al. 20061 . 



Ohkitani fc Kida 1992 . lYeung et al. 1995 



Zhou 19931 ). On the other hand, the and T u b transfers that express the en- 
ergy exchange between the velocity and the magnetic field have a rather different 
behavior. 

4.1 Forced isotropic and homogeneous turbulence 

As mentioned before, forced simulations of MHD turbulence can be attained by 
forcing both fields, or by forcing only the velocity field (in which case magnetic 
fields are sustained by dynamo action; a distinction must be done then between 
the kinematic regime, where the magnetic field has no backreaction on the flow, 
and the turbulent steady state where the magnetic field modifies the flow through 
the Lorentz force) . The mechanically forced case is of more interest as it is closer 
to astrophysical and geophysical configurations, and as it is consistent with the 
constraint of magnetic flux conservation. The first studies of shell-to-shell transfe r 



from simulations in such configuration were presented in lAlexakis et al. (2005bl ). 
In the simulations with a resolution of 256 3 grid points, the velocity field was 
forced with time independent mechanical forcing until a hydrodynamic turbulent 
steady state was reached. Two different forcing functions were studied: one 
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non-helical and one helical. Magnetic Prandtl numbers of unity and smaller 

were considered. Once a hydrodynamic steady state was reached for each forcing 

function, a small magnetic field was introduced and, after the transient kinematic 

dynamo amplification, the system reached a steady state MHD turbulent regime. 

In such a state the transfer functions described in the previous section were 

computed. Typical results are illustrated in Fig. [TJ 

The T uu and T^b transfers were observed to behave in a similar fashion, giving 

direct and local transfer of energy. In Fig. [IJ this is indicated by the negative and 

positive peaks, which show energy is removed by these transfer functions from 

smaller wavenumbers and given to slightly larger wavenumbers. However, for 

T u \, a distinct behavior appeared. The large scale flow injected energy (through 

stretching) directly into the magnetic field at all scales. This manifests itself as 

a peak at the mechanical forcing scale for all receiving shells, and as an extended 

positive plateau (note positive T u b indicates energy given by the velocity field at 

shell Q to magnetic field at shell K). In other words, at a given shell K, the 

magnetic field receives energy from the velocity field in all shells Q < K, and 

gives energy to the velocity field in shells Q > K . This result, reminisc ent of 



the theoretical arguments by 



Batchelor (19501 ) and 



Zel'Dovich et al. (1984 ). was 



interpreted as the sustainment of the magnetic field against Ohmic dissipation 
by dynamo action: to maintain the magnetic field when only the velocity field 
is stirred, a nonzero flux from the velocity field to the magnetic field is required 
at all times. It is worth pointing out here that in the steady state this non-local 
transfer is small compared with the local transfers (approximately 10 — 20% at the 
resolutions studied). When considering Elsasser variables, the transfer functions 
were observed to become more local. 
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The case of random forcing with magnetic Prandtl number of unity was studied 

Carati et al. (20061 ) using 512 3 simulations. The analysis, which used logarith- 



m 



mic binning, confirmed the previous results, showing local transfer in T uu and 
Thb, and non-local coupling between the velocity and the magnetic field. This 
indicates that the phenomenon may be independent of the type of forcing, and 
associated to the stretching process that sustains the magnetic field. The work 
also discussed the possibility of splitting the transfer functions to discriminate 
between forward and backward contributions, which were used to discuss im- 
plications of the shell-to-shell transfers for LES mod els. Similar res ults were 



obtained for forced two dimensional MHD turbulence (IDar et al. 2001 

A different approach was considered by lYousef et al. (20071 ) , who studied the 



steady state of small-scale dynamo action for Pm < 1- Instead of using transfer 
functions, to measure the different components of the energy flux they considered 
the Politano-Pouquet relations © in terms of the velocity and the magnetic field 



5u\\ (|<5u| 2 + |<5b| 2 ))T(^|| (|H 2 + l<5b| 2 ))±2(<5u-<5b (<Su|| t Sb\ 



(22) 



together with 



Chandrasekhar (195ll ) law 



<JuJ) -6(6f<J«||) = --el, 



(23) 



where e is the total energy flux. The authors discriminated between the different 
terms to look how they balanced to give rise to the direct flux. Each of the terms 
in these expressions can indeed be associated to the counterpart in real space of 
the H uu , n^, and TL u b + fluxes in Fourier space. 

The dominant balance was identified between (4/5)eZ and 6 (b^5u^, and they 
concluded that at their available resolution, the local direct cascade of energy was 
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"short-circuited" by the transfer of kinetic energy into magnetic energy. They 

also associated this non-local coupling with the folded str ucture of the small-scale 

magnetic field. Using the shell-to-shell transfer approach, lAlexakis et al. (2007al ) 



further showed that the non-local effects disappear if phases are randomized for 
the two fields, which also make the current sheet and folded structures disappear. 
The non-local effects play a more important role in the kinematic dynamo 



regime ([Mininni et al. 2005al ). as in that case the turbulence is not in a steady 
regime and T u b accounts for all mechanisms that amplify the magnetic field. In 
that case, the T u b transfer has been sh own to be useful to identify and quantify 



scale- by-scale sources of dynamo action (jAlexakis et al. 2007a . iMininni et al. 2005al ). 



4.2 Freely decaying turbulence 

The non-local effects observed in forced turbulence ar e either absent or negligible 
in the freely decaying case. In lDebliquv et al. (20051 ). 512 3 simulations of freely 



decaying MHD turbulence were considered. The T uu and transfers are similar 
to the forced case (see Fig. [2j) and indicate local direct transfer. However, the 
T u b and transfer functions were also observed to be local, with most of the 
transfer between the velocity and the magnetic field taking place between the 
same shell. The remaining transfer (for non-neighboring shells) was observed to 
decay more slowly than in the T uu and T^b functions, but except for this detail 
no other indications of non-locality were reported. 



Similar results were obtained from analysis of solar wind turbulence (jStrumik Macek 201 



Solar wind turbulence is often considered t 
freely decaying wind tunnel turbulence (see 



re MHD equivalent of h ydro dynamic 



Bruno Carbone 20051 for a review) 



From 1996 Ulysses magnetometers time series and using a Markov process ap- 
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proach, 



Strumik &: Macek (200 8b) concluded that the transfer of magnetic to 



magnetic energy was local. Then, using velocity and magnetic field time series 
from ACE spacecraft from 1999 to 2006 and performing the same analysis on the 
remaining transfers they concluded that all transfers were local. 

The differences between the forced and freely decaying cases can be understood 
as in the mechanically forced runs, the velocity field has to continuously supply 
energy to the magnetic field in order to sustain it against Ohmic dissipation. This 
is not necessarily the case for freely decaying runs where both fields are dissipated 
in time. 

4.3 Anisotropic turbulence 

Recently, the shell-to-shell transfers were extended to consider anisotropics when 
an external uniform magnetic field is imposed. This case is of interest as in 
many astrophysical problems a strong large-scale magnetic field is present cre- 
ating small-scale anisotropy. Unlike hydrodynamic turbulence, MHD turbulence 
does not recover isotropy at small scales, and theoretical and numerical results 
indicate anisotropy becomes stronger at smaller scales. 

To study anisotropic transfers, different foldings of th e shells in Fourier spac e 



Alexakis et al. (2007a 



can be implemented. FigOshows the possible options. In 
anisotropic shell-to-shell transfer functions were introduced folding Fourier shells 
in cylinders (associated to wavenumbers k±_ perpendicular to the mean mag- 
netic field) and in planes (associated to parallel wavenumbers fey). Shell-to-shell 
transfers were only considered for the Elsasser variables, but the fluxes were 
reconstructed from these functions to measure the relative contribution of non- 
locality to the total flux. Freely decaying simulations with spatial resolution of 
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256 3 grid points were analyzed, and the amplitude of the imposed magnetic field 

was varied from to 15 (in units of the initial small scale fluctuations). The 
transfer functions of the two Elsasser energies were found local in both paral- 
lel and perpendicular directions, irrespectively of the amplitude of the external 
field. However, interactions between the counterpropagating Alfven waves were 
reported to become non-local. For strong magnetic fields, most of the energy flux 
in the perpendicular direction was found to be due to interactions with modes 
with fcii = (see Fig. In the parallel direction, however, fcii = modes 

cannot transfer energy and most of the interactions were observed to take place 
with modes near k\\ ~ 0. The res ults are in qualitati ve agreement with predic- 



tions from weak turbulenc e theory 



phenomenological models (jAlexakis 20071 ). 



Galt ier et al. 20001 ) and with recent non-local 



Teaca et al. ( 20091 ). 



A different approach to study anisotropic transfers was presented by 
who decomposed the spectral space into rings, studying then transfers along ra- 
dial and angular directions in spectral space (which they termed "ring-to-ring" 
transfers). They considered forced simulations of MHD turbulence with an im- 
posed magnetic field, with a spatial resolution of 512 3 grid points and varying 
the imposed magnetic field between to vlO (in units of the small-scale mag- 
netic field fluctuations). They also observed dominance of energy transfer in the 
direction perpendicular to the uniform magnetic field, and suppression of the 
transfer in the parallel direction. Their approach is useful to understand how 
energy is angularly distributed in spectral space to create anisotropy. Non-local 
effects with the forcing shell were observed in the shell-to-shell transfers, but in 
the angular ring-to-ring transfers were too weak to be noticed. 
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4.4 Magnetic helicity and the inverse cascade 
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Non-local transfers were also reported in investigations of the cascade of mag- 
netic helicity. Magnetic helicity is an ideal invariant in M HD, that is known to 



Brandenburg 200ll. 


1 

Brandenbure & Subramanian 20051. 


Gomez & Mininni 2004. 


Meneeuzzi et al. 1981. 


Pouquet et al. 1976 


). The generation of large scale mag- 



netic fields in galaxies an d other astrophysical bodies is sometimes attributed 



to this inverse cascade. In 



Alexakis et al. (20061 ) magnetically and mechanically 



forced simulations were considered. In both cases, both local and non-local trans- 
fers were observed. At early times, magnetic helicity was observed to cascade 
inversely and locally from the closest neighbor shells, and non-locally from the 
forced shells. When the correlation length became of the size of the box, the 
direct input from the forced scales became dominant, and a local direct trans- 
fer of helicity from large to small scales also developed. This latter effect was 
speculated to be dependent on boundary conditions and therefore non-universal. 
In the mechanically forced ca se the inverse cascade of helicity was associated to 



the large-scale dynamo q-effect ([Br andenburg 2001 



Krause fc Raedler 198(1 



Pouquet et al. 1976 



Branden burg &: Subramanian 20051 . 



Steenbeck et al. 1 966). In that case, 



the mechanical forcing creates equal amounts of magnetic helicity of opposite 
signs at large and small scales. The process can be unde rstood using the concep- 



tual "stretch, twist, and fold" (STF) dynamo mechanism ([Childress fc Gilbert 1995 



Vainshtein &: Zeldovich 1972 ) . Each time a closed magnetic flux tube is twisted 
by the helical velocity field, magnetic helicity is created at large scales, while 
small scale magnetic field lines are twisted in the opposite direction thus creating 



equal amount of magnetic helicity of opposite sign in the small scales. As the 
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STF process is repeated, the large-scale helicity is transferred inversely both lo- 
cally and non- locally (with constant flux), while the small-scale helicity is pushed 
towards smaller scales (see Fig [5]) . This latter process removes magnetic helicity 
from the large scales and allows the magnetic field to "dise ntangle" through re- 



connec tion events, destroying in that way magnetic helicity ([Alexakis et al. 20061 . 



2007bl ). At this moment it is unclear whether this processes should be associated 



to a cascade (i.e., if the p roces s takes place with constant flux), although results 



m 



Alexakis et al. (2007bl ) and iMininni &; Pouquet (20091 ) suggest this may not 



be the case. 

5 Non-local interactions and universality of MHD turbulence 

The above considerations led several authors to consider whether some of the 
usual assumptions in hydrodynamic turbulence hold in the MHD case. From the 
shell-to-shell transfer, the scenario pictured in Fig. [6]seems to arise for the energy: 
interactions between the same fields are mostly local, and interactions between 
the velocity and the magnetic field can have different degrees of non-locality 
depending on whether the turbulence is forced or freely decaying, depending on 
how the velocity and the magnetic fields are maintained against dissipation in 
the forced case, and depending on the presence of an external magnetic field. 
It is unclear for the moment whether the varying degree of non-locality with 
the configuration will converge to a universal solution for very large Reynolds 
numbers. 

Theoretical arguments considering interactions in MHD turbulence also ob- 
tained conflicting results. Using the EDQNM closure, I Pouquet et al. f 19761) re- 



ported non-local interactions which were associated to Alfven waves. In 



Verma f2003l . 



ffcai e Inte ractions in MHD Tu rbulence 



200J) and 
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Verma et al. (20051 ). field-theoretic calculations were used to compute 
the shell-to-shell transfers and it was concluded that they were local except for 
the transfer between the velocity and the magnetic field, which was found to 
be some what non-local. The h elicity transfer was also found to be non-local. 
Recently, 



Aluie &: Evink (2009bl ) gave strict bounds for fluxes in MHD turbu- 
lence under the assumptions that both the velocity and the magnetic energy 
follow power laws in the inertial range between k~ l and k~ 3 . The velocity-to- 
velocity and magnetic-to-magnetic fluxes were found to be local in the limit of 
infinite Reynolds number, and the fluxes coupling velocity and magnetic fields 
were found to be local although counterexamples to their proof as the ones men- 
tioned in Sect. 2 were acknowledged. However, these results shed light into why 
some simulations were found to be more local than others, as mechanisms as 
the small-scale dynamo can be expected to be less relevant in freely decaying 
turbulence in approximate equipartition between the two fields. 

At presently attainable spatial resolutions, other indications of possible non- 



Dmitruk et al. (2003 ). 



universal behavior has been reported in numerical simulations. In[ 
simulations of forced reduced magnetohydro dynamics (RMHD) where presented 
where the energy spectrum changed its power law depending on the timescale of 
the external forcing. Spectra compatible with Kolmogorov, Iroshnikov-Kraichnan, 
weak turbulence theory, or even steeper laws were observed. The RMHD equa- 
tions correspond to an approximation of the MHD equations when a strong exter 



nal magnetic field is imposed. Similar results were reported by 



Mason et al. (2008 ). 



who considered forced MHD with an imposed ma gnetic field. Other numeri- 



cal simulations of forced MHD turbu 



Haugen et al. 200c 



Miiller et al. 2003 



ence (see e.g.. iBeresnvak &: Lazarian 2009 



Miiller Grappin 20051 ) also reported con 
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flicting results. In freely decaying isotropic turbulence, some simulations were 

observe d to develop Iroshnikov-Kraichnan scaling while others Kolmogorov-like 



scaling (jMininni Sz Pouquet 20071 . 



2009, 



Miiller &: Grappin 20051 ). Recently, large 



resolution simulations of freely decaying MHD flows showed that depending on 
the amplitude of the dynam ically consistent large-scale magnetic field, different 



power laws can be realized 



Lee et al. 20091) . Finally, recent studies of spectral 



laws in solar wind data (jPodesta et al. 20071 ) indicate that many of these power 



laws can also be identified in space plasmas. 

Although the main aim of this review is to consider studies of scale inter- 
actions in MHD, in this context it is worth mentioning some of the existing 
phenomenological theories for MHD turbulence. While Iroshnikov and Kraich- 
nan considered small scale fluctuations as isotropic, it is clear now that MHD 



turbulence does not recover isotropy at small scales (jGoldreich &: Sridhar 19951 . 



Milano et al. 2001 



Shebalin et al. 19831 ) and may become even more anisotropic 



as the scales are decreased. To take this into account , a different MHD spec- 
trum has been advocated in lGoldreich &; Sridhar (19951 ). whereby the anisotropy 



of the flow induces a Kolmogorov-like spectrum in the perpendicular direction. 
A balance between linear and non-linear timescales (the Alfven and turnover 
times) is assumed which leads to a "critical balance" of the form /cy-Bo ~ kj_bi. 
Another ani sotropic model based on dynamic alignment of the velocity and mag- 



netic fields (|Boldvrev 20061 ) gives IK-like scaling in the perpendicular direction. 
In this case, the angle between the two fields decreases (and therefore the fields 
become more aligned) with the scale as ~ Z 1//4 . Consideration of this alignment in 
the Politano- Pouquet relations leads to the aforementioned scaling for the energy 
spectrum. Early extensions to flows with sizable cross-correlations can be found in 
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Galtier et al. (200d ) 
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Grappin et al. (19831 ) . Other models have considered tran- 



sitions from Kolmogorov to Iroshnikov-Kraichnan s caling by taking different c om 



binations of the non-l inear and Alfv en timescales (jMatthaeus &; Zhou 19891 ) . or 



consider non-locality (jAlexakis 20071 ) 



Therefore, while the assumptions of locality and of isotropization of the small 
scales common in hydrodynamic turbulence allow for a simpler phenomenologi- 
cal treatment of MHD, th e development of local anisotropies, the variety of time 



scales in the problem (see I Zhou et al. 20041 for a review) , and the different sim- 



ulations showing scaling consistent with different phe nomenological theories, le d 



some authors to discuss some of these assumptions. In lSchekochihin et al. (20081 ) 



non-locality, anisotropy, and non-universality were considered as defining prop- 
erties of MHD turbulence. The authors argued that the small-scale dynamo, a 
fundamental proces s in MHD turbulence, shows clear signatures of non-locality 



Haugen et al. 20031 . 



2004 



Mininni et al. 2005a 



Schekochihin et al. 2002a 



2004, 



2002b). They also argued that anisotropy is intrinsic to MHD, and that non- 
universality manifests itself just from the needed distinction between MHD tur- 
bulence in the presence and in the absence of a strong me an field. Similar 



concerns about universal behavior in MHD were discussed in 



Lee et al. f20Qgh 

for the case of freely decaying turbulence. iBeresnvak &; Lazarian (20oJ ) con- 
sider the lack of a "bottleneck" in MHD (an accumulation of energy at the 
beginning of the viscous range observed in hydrod ynamic turbulence) as evi- 



dence of non-locality (see also 



Graham et al. 20091 ). In some sense, some of 



these discussions can be tracked back to early considerati ons of freely decay- 



ing MHD turbulence and the processes of selective 



Matthaeus fc Montgomery 19801 . 



decay (IKinnev et al. 19951 . 



Mininni et al. 2005b 



Ting et al. 19861 ) and dy- 
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namic alignment 



I I 

Ghosh et al. 198a . brappin et al. 19831 . 
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Mininni et al. 20 05c. 



Pouquet et al. 19861 ). MHD, having three ideal invariants, is known to decay for 



very long t imes into different attractor s depending on t he initial ratio of these 



invariants (jStribling Matthaeus 1991 



Ting et al. 19861 ). Although these solu- 



tions involve final stages of the decay, recent numerical simulations showed that 
these relaxed state s can be realized locally in the flow in very short times scales 



(Mason et al. 2006 



Matthaeus et al. 2008 



Perez k Boldvrev 20091 . 



Servidio et al. 20081 ) 



giving rise to different regimes. 



6 Concluding remarks 

Since the success of the phenomenological theory of Kolmogorov in hydrodynamic 
turbulence, several attempts have been made to apply similar considerations to 
tackle MHD turbulence. The presence of waves and of several time scales, and 
of several ideal invariants limited these approaches giving rise to many possible 
models. Solar wind observations and numerical simulations later showed that 
assumptions like isotropy of the small scales, or equipartition between the fields, 
may not hold in the MHD case. More recently, the increase in computing power 
allowed for some exploration of the parameter space giving rise to conflicting 
results for scaling laws in the energy spectrum. 

The recent introduction of shell-to-shell transfers allowed for detailed studies 
of scale interactions in MHD turbulence and opened the door to discuss another 
hypothesis: that of locality of interactions between scales. The results, at in- 
termediate spatial resolutions and Reynolds numbers, show different degrees of 
non-locality depending on the configuration studied: forced or freely decaying 
turbulence, in presence or in absence of an external magnetic field, etc. Non- 
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local transfers, when observed, involve the coupling between the velocity and the 

magnetic field, or the transfer of magnetic helicity. In the former case, the non- 
local transfers were not larger than 10 — 20% of the total, although they played 
fundamental roles, e.g., sustaining the magnetic field by dynamo action against 
Ohmic dissipation. 

In spite of some conflicting results in the simulations and theory, there is grow- 
ing consensus that MHD turbulence is less local than hydrodynamic turbulence. 
To what extent it is a matter of debate. It is unclear for the moment whether 
these effects will go away for larger Reynolds numbers, or if they stay how much 
impact they will have in the flow dynamics, and under what conditions. How- 
ever, the different degrees of non-locality observed at present resolutions, and the 
existence of non-local processes in MHD (as, e.g., the small scale dynamo), call 
for a discussion about the validity of the hypothesis of locality of interactions, 
and of whether there is only one kind of MHD turbulence or many. This raises 
the question of what is the definition of MHD turbulence in phenomenological 
or theoretical approaches. If only configurations as the ones in the solar wind 
(with an imposed magnetic field) are to be considered, then a universal scaling 
(or several classes of universality) may be possibly identified. However, if pro- 
cesses like the small-scale dynamo, the large-scale dynamo, and inverse cascades 
are to be considered as manifestations of MHD turbulence, non-local interactions 
and non-universal behavior may persist even for very large Reynolds numbers. 
In this context, many of the works reviewed here may have to be revisited in the 
following years, as experiments and increased computing power will allow us to 
explore new regions of the parameter space of MHD turbulence. 



Mininni 

Summary points 

1. Simply applying properties of hydro dynamic turbulence to the MHD case 
may not be possible. In particular, assumptions of scale locality of MHD 
turbulence must be tested in experiments and simulations. 

2. Shell-to-shell transfer functions allow for detailed studies of coupling be- 
tween fields and scales in numerical simulations. The shell-to-shell trans- 
fers can also be associated with physical processes such as as Alfven wave 
interactions, Joule damping, and dynamo action. 

3. The degree of non-locality observed at the presently attainable spatial res- 
olutions depends on the configuration. 

4. Mechanically forced turbulence shows local transfer of magnetic and ki- 
netic energy, but the coupling between the velocity and magnetic field that 
sustains the latter against Ohmic dissipation is non-local. 

5. In freely decaying MHD turbulence, non-local effects seem to be negligible. 

6. Studies of the energy transfer in the presence of an imposed magnetic field 
show most of the transfer takes place in the direction perpendicular to the 
external field, with strong non-local interactions with modes with fey = 0. 

7. The transfer of energy for the Elsasser variables is more local than the 
transfer in terms of the velocity and magnetic fields. 

8. The shell-to-shell transfer of magnetic helicity is more complex, with su- 
perimposed direct and inverse transfers. The inverse transfer has a local 
component, and a non-local one that moves energy from the forced scale 
directly to the largest scales in the system. 
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Figure 1: Left: transfer functions in mechanically forced MHD turbulence, for 
Q = 20. The T uu and functions are local, with a negative peak for K < Q and 
a positive peak for K > Q, which indicate energy is removed by these transfer 
functions from smaller wavenumbers and given to slightly larger wavenumbers. 
The transfer between magnetic to kinetic energy is of smaller amplitude and also 
seems local. Right: The T u f, transfer, for different values of K. This function is 



nondocal, with a strong peak at the forcing scale and with a constant 
plateau that extends up to K w Q. Adapted from lAlexakis et al. (2005bl ). 



positive 



Scale Interactions in MHD Turbulence 



41 



0.1 



-0.1 



Tuu 


Aw 







-5 




K-Q 



Tub 









5 -5 




K-Q 



Figure 2: Left: T uu transfer function in freely decaying MHD turbulence, for 
different shells. The transfer function is similar but has twice the amplitude. 
Right: Tb u transfer function in the same simulation. Note the peak for K — Q = 0, 
indicating most interchange of energy between the velocity and the magnetic 
field takes pla ce between similar sc ales. The shells are logarithmically binned. 
Adapted from iDebliquv et al. (20051 ). 
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Figure 3: Isotropic (spherical) shells in the left, and anisotropic foldings of shells 
in Fourier space. The uniform magnetic field is assumed to be in the z direction. 
Cylindrical and planar shells are shown in the middle, and ring shells are shown in 
the right. Transfer of energy across planes is denoted by T(Qn, Kn), and transfer 
across cylinders is denoted by T(Q±, K±). For ring-to-ring transfers, the notation 
^{k'p} denotes transfer can be measured between K and Q spherical shells, as 
well as between two azimuthal angles a and j3. 
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Figure 4: Left: total energy flux (solid) across cylinders and partial flux associated 
to interactions with modes with /cy = (dashed), with four different values of the 
external magnetic field Bq from to 15 (from top to bottom). Right: same but 
with the total flux and pa rtial flux associated to interactions with modes with 
k\\ = 1 across planes. From lAlexakis et al. (2007al ). 
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Figure 5: Left: the stretch, twist, and fold dynamo mechanism. Each time a 
closed magnetic flux tube is twisted, magnetic helicity of opposite sign is created 
at large and small scales. The folding creates regions where helical magnetic 
fields can reconnect. Right: (a) The helicity spectrum in a simulation with 
(positive) helical mechanical forcing at k = 10. Magnetic helicity is negative at 
larger scales and positive at smaller scales. (b,c,d): The transfer of helicity for 
Q = 2, 10, and 20. The red arrows indicate transfer of negative helicity and the 
green arrows transfer of positive helicity. At large scales (b), negative magnetic 
helicity inversely cascades locally between neighbor shells and non-locally from 
the forced shell and to the small scale shells. At the forced shell (c), the forcing 
injects opposite signs of helicity at large and small scales. At small scales (d), 
positive magnetic helicity has a local direct transfer of helicity, while the small 
scales also remove negative magnetic helicity from the large scales. Note direct 
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Figure 6: Sketch of the several shell-to-shell energy transfers identified in sim- 
ulations of isotropic and homogeneous MHD turbulence. The T uu transfers are 
shown in red, Ty, transfers are shown in blue, and T u b and in green. The 
thickness of the arrows roughly indicates the strength of the transfers. Above: 
mechanically forced simulations. At the shell K, the magnetic field receives en- 
ergy from the velocity field at all larger scales, and gives energy to the velocity 
field at slightly smaller scales. Below: Freely decaying turbulence. The T u f) and 
T{, u transfers only interchange energy between similar scales. In both cases, the 
T uu and transfers are local and give the largest contribution to the flux. 



